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Abstract. We investigate temperature dependent electronic correlation effects in 
the conduction bands of Gadolinium Nitride (GdN) based on the combination of 
many body analysis of the multi-band Kondo lattice model and the first principles 
TB-LMTO bandstructure calculations. The physical properties like the quasi-particle 
density of states (Q-DOS), spectral density (SD) and quasi-particle bandstructure (Q- 
BS) are calculated and discussed. The results can be compared with spin and angle 
resolved inverse photoemission spectroscopy (ARIPS) of the conduction bands of GdN. 
A redshift of 0.34 eV of the lower band edge (T=T'c ^ T=0) is obtained and found in 
close comparison with earlier theoretical prediction and experimental value reported 
in the literature. 
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1. INTRODUCTION 

The rare earth monopnictides and monochalcogenides form a very interesting group of 
materials. They exhibit a rich variety of anomalous physical properties attracting a 
considerable attention for basic and applied research. Among the rare earth alloys, the 
isomorphic compounds like Gd^'^Xp' , where Xp = N,P,As,Sb,Bi and Eu'^^X^^ , where 
Xc = 0,S,Se,Te can be compared since they are isoelectronic in nature. If pure ionic 
bonding is considered then they are expected to be insulators or semiconductors. The 
divalent rare earth monochalcogenides are indeed insulators or semiconductors^ |2] but 
the trivalent rare earth monopnictides, especially the rare earth nitride GdN, show an 
intricate conducting character. For instance, the bandstructure calculations suggest 
GdN to be a semiconductor [3 HI E], half-metal|ni 13 El or a semiconductor in the para- 
magnetic and a semimetal in the ferromagnetic case using quasi-particle self energy 
correct ions [H]. The above mentioned results used different computational techniques. 
Hasegawa et.al, were the first to calculate the self-consistent electronic bandstruc- 
tures of GdN for the paramagnetic case using augmented plane wave (APW) method. 
Their calculations employed one-electron potential instead of local spin density (LSD) 
functional theory. While Lambrecht |1] addressed the problem by estimating the gap cor- 
rections beyond local density approximation (LDA), Ghosh|3] performed self consistent 
spin polarized calculations using full potential linear muffin-tin orbital (FP LMTO) but 
with rigid shifts in the 5d and 4/ states in order to fit the experimental X-ray photoe- 
mission spectroscopy (XPS) and X-ray Bremsstrahlung Isochromat spectroscopy (BIS) 
results. The Self Interaction Corrected local spin density (SIC-LSD) calculations by 
Aerts et.al, ^ and augmented spherical wave (ASW) within LDA and generalized gra- 
dient approximations (GGA) by Eyert^ rendered half- metallic nature to GdN. There 
has also been an interesting investigation on the electronic structure and magnetic prop- 
erties of GdN by Duan et.al, [Zj based on first principles calculations as a function of unit 
cell volume. They found that GdN transforms first from half-metallic to semi-metallic 
and then finally to a semiconductor upon applying stress. These features reveal a strong 
lattice constant dependence on the electronic structure of GdN. However, experimental 
results demonstrated bulk GdN to be a low carrier semimetallicjlOj or GdN thin films 
to be insulating jTTj. 

The magnetic properties of GdN also form a very interesting study and the 
mechanism of magnetism is highly intriguing. There had also been a dispute re- 
garding the magnetic properties jTUj of GdN with earlier reports describing it to 
be an antiferromagnetic|in[ IT2] material while other studies indicated it to be a 
ferro magnet [ini El . But after such controversial discourse, it has been accepted that 
stoichiometric GdN is a ferromagnetic fSJ material. There have been theoretical 
attempts in order to explain the magnetism in GdN using RKKY interaction picture 
by KuznietzlTTl and a model proposed by Kasuva[TH] which is related to the cross pro- 
cesses between the f-d mixing and the f-d exchange interaction. In this paper we do not 
discuss the mechanism of such anomalous behaviour of magnetic properties in GdN but 
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it forms a part of future investigations. As a whole, it generates a huge motivation to 
study such an interesting material and understand the physical phenomenon in similar 
strongly correlated systems. 

Here we combine the many body analysis of the multi-band Kondo lattice model 
alongwith the first principles TB-LMTO spin-polarized bandstructure calculations to 
investigate the temperature dependent electronic correlation effects in the conduction 
bands of GdN. We wish to emphasize that upon using bandstructure calculations we 
are not aiming to find the nature of the electronic ground state of GdN but rather 
take its output as the starting point (input) for our many body theory and compare 
the results for the temperature dependence with that of the experimental ones (redshift 
phenomenon). Very recently a similar study has been reported JHl but using a different 
approach for determining the self-energy |2nj and also for the input in the many body 
part. 

The layout of the paper is as follows. In section |21 we develop the complete multi- 
band Hamiltonian. A self energy ansatz, authentic for low carrier densities, is used 
to solve the model Hamiltonian. In section El we present electronic bandstructure of 
the conduction bands of GdN using tight binding linear muffin-tin orbital (TB LMTO) 
program within LSDA and also obtain the inter-band exchange coupling. In the next 
section |3J the many body model is combined with first principles T=0 bandstructure 
calculations in order to calculate the temperature dependence of some physical proper- 
ties like quasi-particle density of states (Q-DOS), spectral density (SD), quasi-particle 
bandstructure (Q-BS) and red-shift phenomenon in ferromagnetic GdN. In the last sec- 
tion we conclude with the proposition of a spin resolved ARIPS experiment to be 
performed in order to validate our results and thus inspire experimental efforts to be 
carried out in order to study the intriguing temperature dependent properties of GdN. 

2. THEORETICAL MODEL 

The temperature dependent electronic correlation effects in the conduction bands of 
GdN are due to the exchange interaction between the electron in extended {5d) band 
states and the localized (4/) electrons which give rise to localized magnetic moments. 
This situation is covered by the so called ferromagnetic multiband Kondo lattice model 
(KLM), the Hamiltonian of which reads; 



H — Hkin + Hi, 



int 



(1) 



where 



kin / ^ -l^ij '^iaa'^jPcr 
ijafScr 




(2) 
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and 



Hint — ( ^g- Sj cj^^ Cjaa + c\^_^Ciaa) (3) 



The term iffcm denotes the kinetic energy of the conduction band electrons, Sdorbitals in 
case of GdN, with c\^^ and Ciaa being the fermionic creation and annihilation operators, 
respectively, at lattice site Ri. The latin letters (i,j,...) symbolize the crystal lattice 
indices while the band indices are depicted in Greek letters and the spin 



is denoted as a{= t;i)- The multi-band hopping term, T°f is obtained from an 
LDA calculation which incorporates in a realistic manner the influences of all those 
interactions which are not directly accounted for by our model Hamiltonian. It is 
connected by Fourier transformation to the free Bloch energies e°^(k) 

= ^ 5Z ^"''(^) e-*-(^'-^^) (4) 

k 

Hint is an intra-atomic exchange interaction term being further split into two subterms. 
The first describes the Ising type interaction between the z-component of the localized 
and itinerant carrier spins while the other comprises spin exchange processes which are 
responsible for many of the KLM properties. J is the exchange coupling strength which 
we assume to be k-independent and S'f refers to the localized spin at site Ri 

5r = 5r + izA' ; z^ = +l^z^ = -l (5) 

The Hamiltonian in eq. ^ provokes a nontrivial many body problem that cannot be 
solved exactly. The details of the many body analysis along with a model calculation 
and limiting cases are explained elsewhere (see Ref. j2I])- We explain it in brief as 
follows. 

Using the equation of motion method for the double-time retarded Green 
function |22j 

GZSE) = {K^:cI^;))e (6) 

where l,m and are the lattice and band indices respectively, we obtain higher order 
Green functions which prevent the direct solution. Approximations must be considered. 
But, a rather formal solution can be stated as 

G^a{E) = hT[{E + iO+)/ - ?(k) - Sk(T(S)]-^ (7) 

where we exclude the band indices by representing the terms in a generalized matrix 
form on symbolizing a hat over it 

GimAE) = ^ E ^k^(^) e-''^-(^-^-) (8) 

k 

The terms in eq. |71 are explained as follows : J is an identity matrix and e(k) is a 
hopping matrix with the diagonal terms of the matrix exemplifying the intra-band 
hopping and the off-diagonal terms denoting the inter-band hopping. The self energy. 
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Sk(T(E), containing all the influences of the different interactions being of fundamental 
importance, can be understood using site representation : 

im^t, Q,.]-; cL.)) = E K;jaiE)GTmAE) (9) 

Now, we are left with a problem of finding a multi-band self energy ansatz in order 
to compute the Green function matrix and thereby calculate the physical quantities of 
interest like the quasi-particle spectral density (SD) 

A^aiE) = --ImTT{Gkcr{E)) (10) 

and the quasi-particle density of states (Q-DOS) 

k 

According to our many body theoretical analysis|21j, we utilize the multi-band inter- 
polating self energy ansatz (ISA) which is accurately defined in the low carrier density 
regime, for all coupling strengths and satisfying various limiting cases of the model like 
the one of ferromagnetically saturated semiconductor as discussed in the Appendix of 
Ref. 12]- The ansatz is given as : 



(12) 



where 



M^ = z^{S'); a^ = S{S+l)-M^{M^ + l). (13) 
and the bare Green function matrix is defined as : 

G,(E) = ly (14) 

As seen in eq. ^1 we are interested only in the local self energy 

^-(^) = ^E^ka(i?) (15) 

k 

The wave-vector dependence of the self energy is mainly due to the magnon energies 
hw{}i) appearing at finite temperature. But, we neglect a direct Heisenberg exchange 
between the localized spins since we are only interested in the influence of inter-band 
exchange on the conduction band states in order to study the electronic correlations 
and not aimed at calculating the magnetic properties via self consistent calculation of 
localized magnetizations. This can be interpreted as the hx){\<i) —>■ limit. The localized 
magnetization (S^) shall be considered as an external parameter being responsible for 
the induced temperature dependence of the band states. 
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3. Bandstructure calculation 

In order to have the single particle excitation energies i.e the hopping matrices, which act 
as an input in the many body part, and also to find the exchange coupling strength we 
perform the TB-LMTO jlHl bandstructure calculations within LDA. GdN crystallizes 
in a rocksalt structure with an experimental lattice constant of 4.99 A and symmetry 
group Fm-3m. Each Gd^^ ion has twelve nearest and six next nearest Gd-neighbours. 
Figure HI indicates the calculated spin-dependent conduction bands of GdN. Our 




Figure 1. LDA bandstructure of 5d conduction bands of GdN. Full (red) lines are 
spin 1 and dotted (green) ones are spin |. 

evaluation is restricted only to 5 d bands while considering GdN to be a semiconductor 
in accordance with Ref. jSl IH E]- But, even after making such simplifications in our 
calculation, there are typical difficulties due to the LDA which arise in the treatment of 
the strongly localized character of the 4/ levels. In order to circumvent this problem we 
considered the 4/ electrons as core electrons, since our main interest is focused on the 
response of the conduction bands on the magnetic state of the localized moments. For 
our purpose, the 4/ levels appear only as localized spins. Moreover, since we are mainly 
concerned in overall temperature dependent correlation effects, the extreme details of 
the bandstructure are surely not so important. 

In Figure El the LDA- density of states is displayed. A distinct exchange splitting 




Figure 2. Spin-dependent density of states of 5d bands of GdN. Using the center of 
gravity of the bands, the exchange splitting amounts to AE = 1.237 eV. 
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is visible which can be used to fix the interband exchange couphng constant J in eq. El 
Assuming that the LDA treatment of the ferromagnetism is quite compatible with the 
Stoner (mean field) picture, the T=0 splitting amounts to AE=JS. On taking the center 
of gravity of both the bands and alongwith the above assumption we obtain: 

AE = 1.237eV J = 0.35eV (16) 

In the bandstructure calculations, one cannot switch off the inter-band exchange inter- 
action (eq. EI) in the LDA code, but one can speculate from the exact T=0 result for 
empty conduction band|2n] that such an interaction leads only to a rigid shift in the 
t spectrum while the | spectrum is remarkably deformed by correlation effects due to 
spin exchange processes between the delocalized and localized states. So we take from 
the LDA calculation, which holds by definition for T=0, the | part as the single-particle 
input in eq. El Therewith, it is guaranteed that all other interactions which do not ex- 
plicitly come into our Hamiltonian are implicitly taken into account by the LDA single 
particle Hamiltonian. On the other hand, a double counting of any decisive interaction 
is avoided since in the j-spectrum the inter-band exchange (eq. El) only shifts the energy 
zero. Since the methodology used in Ref. [T^ remains the same, it is worthwhile to men- 
tion that the authors have used paramagnetic bandstructure as an input in their many 
body part though it is not guaranteed that the paramagnetic input is free of correlation 
effects. 

In order to get a first impression of the correlation effects we have evaluated our 




Figure 3. The same as in Figure |21 witli additional T=0 results of our combined 
many body and first principles theory for exchange coupling J = 0.35 eV, shown in 
broken lines. 

theory for T=0 K and the results are presented in Figure El The | spectrum is unaf- 
fected by the actual value of J and coincides with the respective LDA curve, after a 
compensation of the unimportant rigid shift (^). The slight deviations, seen in the 
upper part of Figure El are exclusively due to numerical procedures. This shows that our 
theoretical approach for implementing LDA input into the many body model fulfills the 
exact limit (T=0, a = f ) and also circumvents the earlier mentioned double counting 
problem. The lower half, | spectrum, of Figure El demonstrates that correlation effects 
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do appear even at T=0 K. Apart from a band narrowing, they provoke strong deforma- 
tions and shifts with respect to the LDA result. 

4. Model Evaluation 

In this section, we present the numerical results emphasizing the temperature dependent 
electronic correlation effects in the conduction bands of GdN by combining the first 
principles TB LMTO bandstructure calculations within LSDA along with the many 
body theoretical model. 

The correlation effects are examined as the consequence of a test electron created 
( or annihilated ) in an empty (or filled ) band. In our case, we create an electron in an 
empty conduction band with a consideration of GdN as a ferromagnetic semiconductor. 
Such effects are studied by calculating the single-electron Green function (eq. Ej). But, 
apart from electronic sub-system, we also have the magnetic sub-system. The exchange 
coupling between the itinerant electron and localized spins adds up to the correlation 
effects as they produce spin-flip transitions and Ising like interactions in addition to the 
kinetic energy. 

The numerical calculations remain the same as performed in Ref. [21]. For the 
sake of brevity, one can understand it as follows. The single particle output obtained 
from the bandstructure calculations are in the form of Hamiltonian and overlap matrix, 
posing a generalized eigenvalue problem to be solved. In order to employ such matrices 
as an input for the many body calculations, one has to perform a decomposition so as 
to reduce the generalized problem to an eigenvalue problem. Using such an approach, 
one can factorize the Hamiltonian matrices and then obtain new Hamiltonian matrices 
using overlap matrices which can be used directly as hopping integrals, e"(k). This 
factorization and direct use of Hamiltonian matrices circumvent the problem of the 
somewhat ambiguous decomposition of the conduction bandjSEl EZl EHj into single non- 
degenerate subbands which is also used in Ref.|19||. Finally, we use the multi-band self 
energy ansatz (eq. IT^ to compute the Green function (eq. Ej) and therewith calculate 
the spectral densities (eq. HU)) and densities of states feq. ITT|). 

The quasi-particle density of states (Q-DOS) are calculated for J=0.35 eV and 
plotted in Figure^for various values of 4/ magnetizations i.e, different temperature. As 
observed (from the inset of Figure |3)), the lower edge of spin - f Q-DOS performs a shift 
to lower energies upon cooling from T=Tc down to T=0 K which explains the redshift 
effect of the optical absorption edge for an electronic transition. We find a redshift of 
0.34 eV. One should note that a simple mean field treatment with J=0.35 eV and S=3.5 
would result in a redshift of 0.6125 eV. Obviously, the correlation effects drastically 
reduce this value. A similar redshift phenomenon was discussed by Lambrechtjij with 
a calculated value of 0.30 eV and was also reported experimentally having a value of 
0.40 eV by Leuenberger. JH] They observed this effect as a shift of the L2 edge on a 
relative comparison between absoprtion spectra measured on a 2000 A GdN film at T 
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Figure 4. Quasi-particle density of states (Q-DOS) of the 5d conduction bands of 
GdN and as a function of energy for various temperatures (different 4/ magnetizations). 
The outermost curves belong to T=0 K (■^^^=l).On increasing the temperature they 
approach each other. The inset shows on an enlarged scale the temperature shift of 
the lower edge for the case of |. 

> Tc and T < Tc in X-ray magnetic circular dichroism (XMCD) experiment. The L2 
absorption spectra is sensitive to the polarization of t2g states. In Ref.jTHj, the redshift 
of 0.483 eV was obtained using J for the lowest band of the bulk GdN and was found 
to be overestimated. 

Figure El represents the quasi-particle bandstructure for some high-symmetry 

L r XWLKL r XWLK 




L r K Vt L K 



Figure 5. Spin-dependent quasi-particle bandstructure of 5(i conduction bands of 
GdN for different values of magnetization -^g^- 



directions in the first Brillouin zone. The degree of blackening measures the magnitude 
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of the spectral function. As observed in the | spectrum even at T=0 K (^=1), 
parts of the dispersions are washed out showing hfetime effects due to correlation in 
terms of magnon emission and re-absorption with simultaneous spin-flips which can be 
understood as follows. 

At T=0 K and empty bands, the addition of a t electron produces a stable quasi- 
particle since the t electron has no chance to flip its spin with the ferromagnetically 
saturated spin (f|~) sub-system. The imaginary part of the self energy vanishes indicating 
infinite lifetimes. The situation is different for the addition of an | electron as the | 
electron can exchange its spin with the ferromagnetically saturated spin sub-system and 
have finite lifetime. It can emit a magnon and become a t electron provided there exist 
t electron states which can be occupied after spin-flip process. These are the scattering 
states which occupy the same energy regions as the | spectrum. The other possibility 
for the I electron to exchange its spin could be by repeated emission and absorption of 
magnons, thus forming a quasi-particle called magnetic polaron. For strong coupling, 
the above mentioned elementary processes can even lead to correlation- caused splitting 
of the dispersion. This is clearly visible at certain parts in Figure 

For finite temperatures, the spin sub-system is no longer perfectly aligned. There 
are magnons in the system that can be absorbed by the itinerant charge carriers 
provoking quasi-particle damping in the t spectrum too. The spectral weight gets 
redistributed due to spin flip term in the exchange interaction with deformations in the 
density of states. The overall exchange splitting reduces with increasing temperature, 
until in the limit T — Tc {{S^) 0) , the vanishing 4/ magnetization removes the 
induced spin asymmetry. 

In order to have a closer look at the correlations, the spin and k -dependent spectral 
densities are plotted at some high-symmetry points (W,K,X) in the first Brillouin zone 
and for different temperatures. While the Q-DOS refer to angle averaged photoemission 
experiment, the k-dependent spectral densities refer to the angle resolved part. 

As seen in Figure well defined quasi-particle peaks appear with an additional 




Figure 6. Spectral density at X - point. The case of -^^^^ = 1.0 is shown by dotted 
Hnes, while = 0.50 is denoted as broken lines and the full lines represent = 
0.0 

spin split below T^. In accordance with the quasi-particle bandstructure, it is observed 
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in Figure ini that four sharp peaks appear at T=0 K for f spectrum and though already 
strongly damped (at lower energies), the same peak sequence comes out in the | 
spectrum. The exchange splitting collapses for T — Tc and the quasi-particle damping 
increases with increasing temperature. For T=Tc (full lines), the spectrum is washed 
out at higher energies exhibiting strong correlation effects. 

Similar explanations hold for the other high-symmetry points (Figure [7| and 
Figure IHl). There's a strong damping as seen along the entire energy spectrum. And 
some of the peaks (lower energy peaks in Figure Ej) are so strongly damped that they 
may not be visible in the inverse photoemission experiment. Altogether, the 5d spectral 
densities exhibit drastic temperature dependencies, what concerns the positions and the 
widths of quasi-particle peaks. 
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Figure 7. The same as in Figure but for W - point. 
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Figure 8. The same as in Figure|51but for K - point. 



5. SUMMARY and CONCLUSION 

We investigated the temperature dependent electronic correlation effects on the 5d con- 
duction bands of GdN. In this respect, we started with multi-band Kondo lattice Hamil- 
tonian which describes an exchange interaction between the localized moments and itin- 
erant conduction electrons. Using the equation of motion approach in terms of Green 
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function, we tried to evaluate the multi-band KLM Hamiltonian. In order to find the 
solution, i.e. Green function, there was a need to find a self energy ansatz which incor- 
porated all the interactions of the system. We used an Interpolating Self energy Ansatz 
(ISA)|23] for the multi-band j25 case which was found to be reliable for low carrier 
densities and all coupling strengths and which fulfilled important limiting cases of the 
model. 

We further performed bandstructure calculations using TB LMTO program in or- 
der to obtain hopping matrices which served as an input in the many body analysis. 
Though there has been a wide discrepancy regarding the nature of the electronic ground 
state of GdN, we considered it to be a semiconductor and thereby concentrating only 
on the empty 5 c? bands and the correlation effects due to their interaction with 4/ lo- 
calized moments. From the bandstructure calculations, we also obtained the inter-band 
exchange coupling strength as a consequence of exchange based splitting of the 5d con- 
duction bands. 

The first principles band structure calculations were combined with the many body 
theory in order to calculate the Green function and therewith some of the physical prop- 
erties of interest in GdN like quasi-particle spectral densities (Q-SD) and densities of 
states (Q-DOS) are estimated for different values of magnetization and exchange cou- 
pling strength. A redhift of 0.34 eV for exchange coupling J=0.35 eV was calculated 
and found to be in close comparison with other theoretical prediction and experimen- 
tal value reported in the literature. Strong temperature dependent correlation effects 
were observed over the whole energy spectrum. Our results can be compared with spin- 
dependent ARIPS of the conduction bands of GdN and thus we would like to instigate 
experimental efforts in order to study this interesting compound. 
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